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Background Of The Invention 

1 . Field of the Invention. 

This invention relates to the estimation, derivation and use of a signal-to noise- 
ratio (SNR) - related parameter, such as Eg/No, for symbols received over a 
communication channel. 

2. Related Art. 

In wireless or wireline communications systems, it is often necessary or 
desirable to estimate a signal-to-noise ratio - related parameter for symbols received 
over a communications channel. One such parameter is Es/No, where Es is the 
received energy/symbol, and No is the noise power spectral density. Another is 
Eb/No, where Eb is the energy^it, and No is again the noise power spectral density. 
The parameter Eb/Nq is related to Es/Noas follows: 
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where Eb and No are as defined previously, and Rbitspersymb is the number of source bits 
(pre- error correction coding) delivered by each channel symbol. 

Both Eb/No and Es/No bear a relationship with SNR and thus are properly 
characterized as SNR related parameters. For example, the parameter Eb/Nq bears the 
following relationship to SNR: 

N, n{r,) 

where Eb and No are as defined previously, S is the average signal power, N (=NoW) 
is the average noise power, S/N is the signal to noise ratio (SNR), W is the system 
bandwidth, and Rr is the bit transmission rate. 

Moreover, the variance, cr^ , of the additive noise borne by the received 
symbols is related to the noise spectral density by the relationship 

Before the symbols are input to a trelhs decoder, for example, it may be 
necessary to estimate an SNR-related parameter such as Es/Nq. The reason is that the 
computation of branch and state metrics performed during the decoding process may 
may be weighted and/or quantized differently depending on the relative value of the 
noise energy associated with the symbols being decoded. What's more, these signal- 
to-noise ratio estimates may be used to assist a user in pointing an antenna to ascertain 
which pointing direction leads to most effective reception. 
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For additional information on trellis decoders, including maximum a posteriori 
(MAP) decoders, log-MAP decoders, Max-Log-Map decoders, Viterbi decoders, and 
Soft Output Viterbi (SOVA) decoders, please see A, Viterbi, "An Intuitive 
Justification and a Simplified Implementation of the MAP Decoder for Convolutional 
Codes," IEEE Journal On Selected Areas In Communications, Vol. 16, No. 2, Feb. 
1998, pp. 260-264; S. Benedetto et al., "A Soft-Input Soft-Output Maximum A 
Posteriori (MAP) Module to Decode Parallel and Serial Concatenated Codes," TDA 
Progress Report 42-127, Nov. 15, 1996, pp. 1-20; D. Divsalar et al., "Turbo Trellis 
Coded Modulation with Iterative Decoding for Mobile Satellite Communications," 
Proc. hit. Mobile Satellite Con£, June 1997; "A Comparison of Optimal and Sub- 
Optimal MAP Decoding Algorithms Operating in the Log Domain," Proc. IC '95, 
Seattle, WA 1995, pp. 1009-1013; C. Berrou et al., "Near Shannon Limit Error- 
Correcting Coding And Decoding: Turbo-Codes," Proc. IEEE Int. Conf Commun., 
Geneva, Switzerland, 1993, pp. 1064-1070; L.R. Bahl et al., "Optimal Decoding of 
Linear Codes For Minimizing Symbol Error Rate," IEEE Trans. Inform. Theory, Vol. 
IT-20, pp. 284-287, 1974, each of which is incorporated by reference herein as 
through set forth in full. 

One approach for estimating total power involves analyzing Automatic Gain 
Control (AGC) settings at the receiver. However, this approach is cumbersome 
because an estimate of total received power does not easily translate into an estimate 
of an signal-to-noise ratio related parameter such as Es/Nq. For one, it may also be 
impractical since extremely precise calibration of the AGC and other receive chain 
gain elements over vohage and temperature variations may not be easy. Moreover, 
knowledge of the noise figure of these same elements may not be characterized, 
especially over temperature variations. Another approach for estimating total power 
is to estimate the noise power, Np, by (1) computing the sample variance of 



U.S. Express Mail No.: EL 711909771 US 
Filing Date: May 30, 2001 



PATENT 
HOW01827.0050.00US00 
00CXT0374D 



yjl^ , (2) estimating the 'signal power + 



noise power', Sp + Np by computing 



the mean of + g% and forming the appropriate ratios and subtractions to derive the 
signal to noise ratio Sp/Np, and then (3) converting Sp/Np to the Es/Nq using the 
relations Sp = Eg and No/2=Np. Here I and Q are the in-phase (I) and quadrature phase 
(Q) components of quadrature symbols received at the receiver. However, this 
computational process is involved. 



A system of and method for estimating an signal-to-noise-ratio (SNR) - related 
parameter, including but not hmited to Es/No, for symbols received over a 
communications channel is provided by the invention. The symbols are analyzed to 
determine if they fall within one or more predetermined collection areas defined in 
relation to a constellation of possible transmitted symbols. In one implementation, 
AGC has been performed on the symbols, so that the average signal power is known. 
The number of received symbols within a group which fall within the one or more 
collection areas is counted. The resulting count may be filtered and/or averaged over 
multiple symbol groups. A predetermined lookup table may be provided to correlate 
values of the count with values of the SNR - related parameter. The SNR - related 
parameter in the table may be expressed in dB in order to increase the dynamic range 
of the table. 

hi one application, a count is empirically determined, and correlated with a 
value of the corresponding SNR - related parameter. This value, if expressed in dB, is 
converted into linear form to form a scale factor which may be used to scale the 
received symbols before they are input to a trellis decoder. The scaled symbols may 
then be quantized. A uniform quantization procedure may be performed. The 
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quantization procedure may be characterized by a quantization delta which has been 
optimized around a nominal value of the SNR - related parameter. 

The count may be updated or refreshed over time. Consequently, the scaUng 
factor may change over time as the count changes. 

5 Other systems, methods, features and advantages of the invention will be or 

will become apparent to one with skill in the art upon examination of the following 
figures and detailed description. It is intended that all such additional systems, 
methods, features and advantages be included within this description, be within the 
scope of the invention, and be protected by the accompanying claims. 

1 0 Brief Description of the Figures 

The invention can be better understood with reference to the following figures. 
The components in the figures are not necessarily to scale, emphasis instead being 
placed upon illustrating the principles of the invention. Moreover, in the figures, like 
15 reference numerals designate corresponding parts throughout the different views. 

Figure 1 illustrates an example of a trelUs representation for a log-MAP 
decoder. 

Figure 2 is a block diagram of a system for scaling received symbols 
according to the invention. 

20 Figure 3 is a flowchart of a method of scaling received symbols according to 

one application of the invention. 
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Figure 4 illustrates an example point on a BER vs. Es/No curve for optimizing 
the quantization delta in the system of Figure 2. 

Figure 5 illustrates one example of quantization levels in the system of Figure 

2. 

5 Figures 6A-6B are flowcharts of a method of estimating an SNR - related 

parameter according to the invention. 

Figure 7 A illustrates an example of a collection area centered within an 8-PSK 
symbol constellation. 

Figures 7B - 7C illustrate examples of uni-dimensional collection areas within 
10 an 8-PSK symbol constellation. 

Figure 7D illustrates an example of off-center collection areas within an 8- 
PSK symbol constellation. 

Figure 7E illustrates an example of a collection area centered within a QPSK 
symbol constellation. 

15 Figure 8 illustrates the process of measuring conditional probabilities in an 

example collection area for an 8-PSK symbol constellation. 
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Detailed Description 

I. Example Application 

In this section, for illustrative purposes only and without limiting effect, an 
5 example application of the invention will be described. However, it should be 
appreciated that many other apphcations of the invention are possible. In particular, 
although the application described here relates to the use of the invention in the 
context of a log-MAP decoder, it should be appreciated that the invention may be 
used in other contexts, and with other trellis decoders besides log-MAP including, 
10 without limitation, soft output decoders, MAP decoders, Viterbi decoders, SOVA 
decoders, and MAX-log-MAP decoders. 

A portion 100 of a trellis representation for a log-MAP decoder is illustrated in 
Figure 1 . The log-MAP decoder may be the inner decoder in a serial concatenated 
convolutional code decoder. The trellis representation illustrates the states of a 
15 corresponding log-MAP encoder and transitions between the states over discrete 
points in time. In Figure 1, numeral 102 identifies the possible states of the encoder at 
time k, and numeral 104 identifies the possible states of the encoder at time k+1. 

In the particular example illustrated in Figure 1, for illustrative purposes only 
and without any hmiting effect, there are three possible states of the encoder 
20 illustrated at time k, mk\ mk^, and mk^, and three possible states of the encoder 
illustrated at time k+1, mk+i\ mk+i^, and m^+i^. The branches between the states 
represent possible hypotheses about the symbol that was transmitted at time k. 
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During the decoding process, state metrics are computed for each of the states 
in the trellis representation, and branch metrics are computed for each of the 
transitions between the states in the trellis representation. 

State metrics are computed recursively. Both forward state metrics and 
5 reverse state metrics are computed. The forward state metrics are metrics which, for 
states at time k, are computed recursively based on the state metrics for time k-1. The 
reverse state metrics are metrics which, for states at time k, are computed recursively 
based on the state metrics for time k+1. In Figure 1, the forward state metrics for 
states mk^ mk^, and mk^ are respectively identified as ak\ ak^, and ak^, and the reverse 
10 state metrics for mk+i\ mk+i^, and m^+i^ are respectively identified as (3k+/, Pk+i^, and 
Pk+i'. 

A branch metric for a transition between states m and m' at time k is referred 
to in Figure 1 as yk*"™ • A branch metric at time k represents the conditional 
probability of the received symbol at time k, rk, given the transmission of a 

15 hypothesized symbol Sk, at time k. The transitions illustrated in Figure 1 comprise the 
following: 1.) a transition between state 1 at time k and state 1 at time k+1, having a 
branch metric yk^^ representative of the conditional probability Pr(rk I sg); 2.) a 
transition between state 1 at time k and state 2 at time k+1, having a branch metric 
yk^^ representative of the conditional probability Pr(rk I se); 3.) a transition between 

20 state 2 at time k and state 1 at time k+1, having a branch metric yk^' representative of 
the conditional probability Pr(rk I S3); 4.) a transition between state 2 at time k and 
state 2 at time k+1, having a branch metric yk^^ representative of the conditional 
probability Pr(rk I S4); 5.) a transition between state 2 at time k and state 3 at time k+1, 
having a branch metric yk^'^ representative of the conditional probability Pr(rk | S3); and 

25 6.) a transition between state 3 at time k and state 3 at time k+1, having a branch 
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metric yk^^ representative of the conditional probability Pr(rk I S5). For illustrative 
purposes only, these transitions are assumed to correspond respectively with the 
following hypothesized symbols: Sg, Se, S3, S4, S3, and S5. 

Several examples of the recursive calculation of the forward state metrics can 
5 be described in relation to Figure 1 . One example is the calculation of the state metric 
aic+i^ This calculation can be expressed using the following equation: 

al,=MAX*{al+rl\al+rl') (1) 
where the function MAX (A,B) is defined to be: 

MAX* {a, B) = MAX{A, B) + hi(l + exp(- 1^ - 5|)) (2) 
10 Another example is the calculation of the state metric ak+i^. This calculation can be 
expressed using the following equation: 

al, = MAX* [al + yf , al + yf ) (3) 
A third example is the calculation of the state metric ak+i^. This calculation can be 
expressed using the following equation: 
15 al, ^MAX*[al + rl\al + ) (4) 

Several examples of the recursive calculation of the reverse state metrics can 
also be described in relation to Figure 1 . One example is the calculation of the state 
metric p]c^ This calculation can be expressed using the following equation: 

/3l = MAX* [pl, + yl' , /3l, + yf ) (5) 
20 Another example is the calculation of the state metric This calculation can be 
expressed using the following equation: 

PI =MAX*{pl,^yl\Pl, -^yT^PL ^yL) (6) 
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A joint probability is the probability of transitioning from a particular state m 
to a state m' at time k over a particular branch. In a log-MAP decoder, due to the 
expression of the state and branch metrics in the natural log domain, which allows 
linear domain multiplication to be performed through addition in the natural log 
domain, the joint probability of transitioning from a state m to a state m' at time k, 
A,k™™ , may be expressed ak™ + yk™™' + Pk+i™'. 

In a log-MAP decoder, these joint probabilities are determined once the state 
and branch metrics have been. Once these joint probabilities have been determined, 
for each possible hypothesized symbol at time k, Si, a log-hkelihood fimction, LLk(Si), 
may be determined in accordance with the following equation: 

LL, (^, ) = M4X;,^,, ,,,, „^ ^ {XT' ) - MA^r;,^,, (XT' ) (7) 

where the ftmction MAX*(A,B) is as set forth above in expression (2). Note that this 
expression takes account of the possibility that more than one transition at time k may 
imply release of the same symbol. This occurrence is illustrated in Figure 1, where 
the branch from state 2 to state 1, and that from state 2 to state 3, both imply the 
release of S3. 

Once the log- likelihood functions have been determined, the decoder then 
selects the hypothesis Sj' at time k which has the largest LL value. This condition 
may be expressed as follows: 

LL^ [s. ' ) - MAX^. {LL^ (s^ )) (8) 

In the case in which the possible symbols s conform to an 8-PSK symbol 
constellation, as explained in co-pending U.S. Patent Apphcation Serial No. To Be 
Determined, Howrey Dkt. No. 01827.0049.00US00, which is hereby fully 
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incorporated by reference herein as though set forth in full, a branch metric for a 
branch in the trellis corresponding to a received symbol r and a hypothesized symbol s 
maybe expressed as follows: 



5 where * indicates the complex conjugate operation, r and s are both complex symbols, 
and is the variance of noise introduced due to transmission through a 
communications channel. 

If r and s are quadrature symbols, r can be expressed as ri + jrg, where ri is the 
in-phase component of r, and rq is the quadrature component of r, and s can be 
10 expressed as Si + jsg, where Si is the in-phase component of s, and Sq is the quadrature 
component of s. Expressing r and s in these terms, the expression rs' + r*s can be 
expanded to rjsi - jriSg + jrqSi + TqSq + rjSi +jriSQ -jrgSi + rqSg, which simplifies to 2riSi 
+ 2rQSQ. Substituting this expression for the numerator of expression (9) yields the 
following for the branch metric: 



In order to ensure proper calculation of the state metrics reflected in 
expressions (1), (3)-(6) above, the denominator cr^ must be included in the branch 
metrics utilized in these calculations. That is because these calculations involve use 
of the MAX* operation, shown above in expression (3), and the MAX* operation, 
20 being non-hnear, requires that the denominator be included so that a relative 
comparison of alternate branches can properly be performed. 
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To take account of the denominator a^, a preferred approach is to scale each of 
the received symbols r, or the components thereof, by a suitable scaling factor, prior 
to using these sjonbols in the computation of the branch metrics. The scaled symbols, 
or the components thereof, would then be combined with the hypothesized symbols s 
5 in the manner prescribed by expressions (9) or (10) above to compute the branch 
metrics. An estimate of the SNR - related parameter such as Es/No would be a 
suitable scaling factor since it acts as a proxy for the parameter \J<3^, assuming that 
the AGC has scaled the signal such that the average signal power Es = 1. The 
following expression defines the relationship between 1/cy^ and Es/No (expressed in 



An alternate approach for taking account of the parameter cr^ is to avoid 
scaling the received symbols, but to swap amongst alternate lookup tables for 
performing the MAX operation, where each lookup table reflects a different value of 
15 a^. However, this approach is not preferred, since most log-MAP decoders require 
performing a large number of MAX* operations in parallel, and to avoid contention, 
most if not all of the MAX* operations will require their own lookup tables. 
Consequently, the amount of swapping that would have to be performed as 
changes might be prohibitive. 

20 Figure 2 illustrates a system 200 for scaling received symbols before inputting 

the scaled symbols to a log-MAP decoder. As illustrated, symbols 202 are received 
over a communications channel and input to the system. The symbols are input to 
multiplier 206 which multiplies the symbols by scaling factor 204. The scaling factor 
is derived from an estimate of an SNR - related parameter such as Es/No. The scaled 



10 dB): 




(11) 
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symbols are then input to quantizer 208, which quantizes the scaled symbols. The 
quantized, scaled symbols 210 are output by the system. They may then be input to a 
log-MAP decoder where they are used in the computation of branch metrics. 

A method of operation of the system of Figure 2 is illustrated in Figure 3. As 
illustrated, the method begins with step 302, where symbols received over a 
commimications channel are scaled using a scahng factor derived from an estimate of 
an SNR - related parameter such as Es/Nq. Step 302 is followed by step 304, where 
the scaled symbols are quantized. 

The scaled symbols may be quantized using a uniform quantizer in which the 
delta A between successive levels in the quantization scale is fixed, but is optimized 
around a certain point on a Bit Error Rate (BER) v. Es/Nq curve. Figure 4 illustrates 
a family of BER vs. Es/No curves 402, 404, and 406, each corresponding to a certain 
Forward Error Correction (FEC) code encoding rate. Numeral 410 identifies a 
maximum BER level which must be maintained in operation. Numeral 408 identifies 
an example BER vs. Es/No point around which the quantization delta A may be 
optimized. Typically, this point will be at about an intersection with the maximum 
BER curve 410. An operating point with an Es/No much beyond then is unnecessary 
since the maximum BER is all that is required to be achieved for reliable 
communication. An operating point with an Es/No less than this is not acceptable 
because it implies a BER greater than the maximum BER which is allowed. 

Once determined, the quantization delta may be fixed regardless of the 
encoder rate. Therefore, the quantization delta, optimized around the point 408 on 
curve 404 in Figure 4, may be used even if a different encoder rate than that 
corresponding to curve 404 is used. 
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In one example, a total of 4 bits is allocated to representing I and Q values of 
received symbols. The result is that each of the 1 and Q values can take on one of 16 
values. The situation is depicted in Figure 5. 

The encoder rate in this example is 2/3, the modulation type 8-PSK, and the 
5 maximum allowed BER correlates with a value of Es/No equal to 8 dB. Optimization 
was performed at a value of Es/No equal to 8 dB. The optimal step size delta, 5, was 
determined using a search algorithm which sought to determine the value of delta that 
maximized channel capacity at a value of Es/No equal to 8 dB. The resulting step size 
delta was determined to be 0.158. 

10 To account for scaling, an extra bit (5 bits total) was allocated to the 

representation of the scaled I and Q values. The scale factor, 1/a^, for a value of 
Es/No equal to 8 dB, was determined using expression (11) to be 12.6191. The 
resulting quantization delta, A, was thus determined to be 0.158«12.6191 = 1.9938. 

II. Preferred Embodiments 

15 An embodiment of a method 600 for estimating an SNR - related parameter, 

including but not limited to Es/Nq, for symbols received over a communications 
channel is illustrated in Figures 6A-6B. The method begins in step 602, where a 
symbol r is received after transmission over a commimications channel. It is assumed 
that the signaling constellation has been 'centered' by AGC, so that the average power 

20 of the signal is known, and set to a fixed value, regardless of signal power transmitted 
or received (pre-AGC). This allows the measurement of what amounts to a 'noise 
spread' to be invariant of the received power. In addition, it is assumed that the 
symbol constellation is stationary (i.e., not rotating or jittering), which implies that the 
demodulator's symbol-time, phase and frequency tracking loops are locked, and are 
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successfully tracking the signal as it continuously varies due to channel variations, 
etc. 

Step 604 follows step 602. Li step 604, the location of the symbol in terms of 
a predetermined coordinate system is determined. Step 606 follows step 604. In step 
606, the location of the symbol is analyzed to determine if it falls within one or more 
predetermined collection areas. 

The one or more predetermined collection areas of step 606 should be defined 
based on the symbol constellation being used. Figure 7A illustrates an example 
collection area 702 for an 8-PSK symbol constellation 700. In this particular 
example, the collection area is centered on the I-Q plane. It is also two-dimensional, 
i.e., defined in terms of both I and Q, e.g., |I| < 2 and |Q| < 2. 

A received symbol r in this example is assumed to be a quadrature symbol 
having in-phase (I) and quadrature phase (Q) components. Step 604 is performed by 
first locating the symbol on the I-Q plane, and determining if the location of the 
symbol falls within the collection area 702. 

Note that the collection area need not be two-dimensional, i.e., defined in 
terms of both I and Q, but can also be uni-dimensional, i.e., defined in terms of just I 
or just Q. Moreover, the collection area need not be centered on the I-Q plane, but 
can be off center. Furthermore, there can be multiple collection areas. Figure 7B 
illustrates an example involving an 8-PSK symbol constellation where the collection 
area 702 is imi-dimensional, and defined in terms of just I, e.g., [I| < 2. Figure 7C 
illustrates an example involving an 8-PSK symbol constellation where the collection 
area 702 is uni-dimensional, and defined in terms of just Q, e.g., |Q| < 2. Figure 7D 
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illustrates an example where two collection areas, identified with numerals 706 and 
708, are defined, each of which are mirror images of the other. 

Also, a collection area is not limited to being defined for an 8-PSK symbol 
constellation, but can be defined for any symbol constellation including, without 
limitation, QPSK, 16-QAM, and 64-QAM. Figure 7E illustrates an example where 
collection area 722 is defined for QPSK symbol constellation 720. 

Turning back to Figure 6A, step 606 is followed by step 608, where the 
location of the received symbol on the predetermined coordinate system is analyzed to 
determine if it falls within the one or more predetermined collection areas which have 
been defined. If so, step 610 is performed. In this step, a value Count, which 
indicates the number of symbols within a group defined by Window_Size which fall 
within one of the collection areas, is incremented. If not, step 610 is bypasssed. 

Step 612 is performed next. In step 612, a value Index is incremented. Index 
is a running total of the number of symbols which have been received in the current 
group. Step 614 follows step 612. In step 614, the value Index is analyzed to 
determine if it exceeds a predetermined value Window_Size. If not, the method loops 
back to step 602. If so, the method continues with optional step 616. 

At the inception of optional step 616, the value Count indicates the number of 
received symbols which have fallen into one of the collection areas out of a total 
number of symbols defined by the value Window_Size. In optional step 616, the 
value of Count is averaged or filtered to remove short term fluctuations. In one 
embodiment, the value of Count is averaged with previous values of Count (stored in 
optional step 620) to form a moving average. In one implementation, the value of 
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Count is averaged with previous values of Count simply by adding the values 
together. 

Step 616 is followed by step 618. In step 618, the value of Count, after any 
averaging or filtering from step 616, is converted into an estimate of an SNR - related 
5 parameter. In one embodiment, the SNR - related parameter which is estimated is 
Es/No- In one implementation of this embodiment, step 618 is performed by 
accessing one or more lookup tables which associate a value of Count with an 
estimate of Es/Nq. (Values within the lookup tables may be tabulated via computer 
simulation.) Once determined, the value of the SNR - related parameter is then 
10 output. 

The association between count values and values of the SNR — related 
parameter is highly dependent on the specific circumstances, but the general principle 
involved is as follows. First, it is assumed that the noise introduced through 
transmission of a symbol through a communications channel is additive, and has a 

15 Gaussian probabihty distribution function (pdf) with zero mean. Second, that since 
the noise has zero mean, and follows a Gaussian distribution, the variance 
completely specifies the noise distribution. Third, the received signal has variance 
equal to the noise variance, and mean equal to the normalized transmitted symbol. 
Fourth, the conditional probablity that a received signal falls within a particular 

20 collection area assuming a certain transmitted signal, can be determined by 
subtracting the assumed transmitted signal from the range of received signals (to form 
a range for the noise residual signals), and then evaluating the probability that the 
noise distribution (with variance a^) lies within that area. Fifth, it is assumed that all 
symbols are, on average, transmitted with equal probability. Sixth, the total 

25 probability that a received signal falls in a particular collection area is the sum of the 

aforesaid conditional probabilities, each multiplied by the assumed symbol 
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probability— which is the same for every assumed symbol. Seventh, a mapping can be 
constructed between this total probabiHty and the noise variance, so that the aforesaid 
probability calculation can be reversed, thereby enabling the noise variance to be 
inferred from the measured total probability — via table lookup, for example. Eighth, 
5 a coimt of the relative frequency (i.e., the number of counts falling within a collection 
area per total number of evaluated received symbols) is a proxy estimate of the actual 
probability. Ninth, if the signal level is normalized, a l/cr^ transformation of a 
estimate yields an SNR-related parameter estimate. Therefore, a value of the SNR - 
related parameter can be determined from relative counts of symbol data falling with 
10 a given collection area. 

In the 8-PSK case, the process of determining conditional probabilities can be 
further explained with reference to Figure 8, which illustrates a uni-dimensional 
collection area 724 (defined in terms of the I axis) superimposed on an 8-PSK symbol 
constellation, and the orthogonal projections (along the I axis) of the conditional pdfs 
1 5 for each of the eight symbols in the constellation. The probability of the received 
signal ri falling into the collection area is the sum of the areas of each of the pdfs 
within the collection area weighted by 1/8. 

During a calibration phase, a table associating the parameter with the 
probability of the received signal falhng within the collection area can be determined. 
20 Then, during operation, a count of the relative frequency of symbols falling within the 
collection area can be empirically determined, and used to lookup the associated value 
of G-^ (or equivalently 1/a^) in the lookup table. 

Step 618 is followed by optional step 620. In optional step 620, the value of 
Count, prior to the averaging and/or filtering of step 616, is stored for use in forming a 
25 moving average in subsequent invocations of step 616. 
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Step 620 is followed by step 622. In step 622, Count and Index are both reset 
to zero. The method then loops back to step 602. 

In one implementation example, a collection area known as a bin is centered 
on the I-Q plane. Symbols within a predetermined count window are analyzed to 
5 determine the number that fall within the bin. The resulting count value is averaged 
over a predetermined averaging window to result in an averaged count value which is 
output. 



The following pseudo-code illustrates the process of outputting an averaged 
count value in this implementation example: 

10 Symbol_Count++; 

if (I <= bin size & I >= -bin_size & Q <= bin size & Q >=bin_size) 
Count++; 

if (Symbol_Count >= Count Window) 
Smn = Sum + Count; 
15 Count = 0; 

SymbolCount = 0; 
Num_Windows++; 
if (Num_Windows >= Averaging_Window) 
Output (Sum); 
20 Sum = 0; 

Num_Windows = 0; 



In this pseudo-code, the size of the bin is defined by bin_size, which specifies 
the extent to which the bin extends from the origin along both the I or Q axes. The 
size of the count window in terms of symbols is defined by Count Window. The size 
25 of the averaging window in terms of count windows is defined by 
Averaging_Window. The count of the number of symbols falling within the bin since 
the last count window ended is defined by Count. The count of the number of 
symbols received since the last count window ended is defined by Symbol_Count. 
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The count of the number of count windows since the last averaging window ended is 
defined by the parameter Num_Windows. The running total of the number of 
symbols falling within the bin since the end of the last averaging window is defined 
by Sum. 

5 The value of Sum is output firom the foregoing process. At this point, Sum is 

then associated with an estimate of Es/Nq. In one implementation example, this is 
achieved by maintaining two lookup tables, Tablel and Table2, where an entry in 
Tablel, Tablel [k], is the probability of a noisy symbol falling within a bin of the size 
bin_size at the value of Es/No indicated by a corresponding entry, Table2[k], in 

10 Table2. Note that Tablel is assumed to be monotonically decreasing, while Table2 is 
assumed to be monotonically increasing. (That follows fi"om the fact that, as Es/No 
increases, the spread of the noise distribution decreases, and, so, if the bin collects 
noisy received data in an area lying between noise-fi-ee symbol locations, the 
probability of a symbol falling in the bin decreases.) Also, the values of Eg/No in 

15 Table2 are expressed in terms of dB; and Table2 depends on the symbol constellation, 
with different instantiations of the table being appropriate for 8-PSK and QPSK, and 
the bin size. 

To perform this association, entries in Tablel are first converted to count 
values by multiplying them by a factor A, which is the product of Coimt_Window and 
20 Averaging_Window. The two count values derived fi-om Tablel which boimd the 
value Sum are then determined. These two entries will bear the following relationship 
to Sum: 

A • Table\[k] > Sum >A» Tablel[k + 1] (12) 
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The two corresponding entries in Table2, Table2[k] and Table2[k+1], are then 
retrieved. The estimate of Es/Nq corresponding to Sum is then interpolated from these 
two values as follows: 



if (Sum A«Tablel [Number_of_rows - 1]) 
row = Number_of_rows - 1; 
bO = A«Tablel[row-l]; 
bl = A»Tablel[row]; 

1 5 ESNO = Table2[row- 1 ] + ((bO - Sum)/(bO - b 1 ))»(Table2[row] - Table2[row- 1 ]) ; 

Note that the value ESNO in this implementation example is refreshed every 
Count Window symbols. 

As an alternative to the foregoing, the values in Table 1 may already in the 
form of count values. In this alternative, the two entries in Table 1 which bound Sum, 
20 Table l[k] and Table l[k+l], are determined. Then, the two corresponding entries in 
Table2, Table2[k] and Table2[k+1], are retrieved, and the value ESNO determined as 
follows: 



5 



Pseudo-code illustrating these operations is as follows: 



10 



for (n = 0; n < Number_of_rows; n++) 

if (Sum >= A»Tablel[n+l] & Sum <= A«Tablel[n]) 

row = n+1; 
if (Sum >= A.Tablel[0]) 

row = 0; 



ESNO = Table2[k] + 



Table\[k] - Sum 



• {Table2[k + 1] - Table2[k]) (14) 



{Tablel[k]-Tablel[k + l]) 



The value ESNO may then be converted to a scale factor, SCALE, for use in 

25 the scaling application described at the outset. The conversion is necessary since 
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ESNO is expressed in dB, and SCALE needs to be in linear terms. In one 
implementation example, this is performed through a third lookup table, TableS, 
which maps values of ESNO to values of SCALE. This table implements the 
following function; 



In one implementation example, the resulting value of SCALE is a 12 bit 
number in the range [0.5, 2.5]. The incoming I and Q components of the received 
symbols are each 10 bits. The resulting scaled symbol (post-scaling multipHcation) is 
21 bits for each component. The scaled symbols are then input to a 5-bit quantizer 
which implements a uniform quantization scheme, where the quantization delta. A, is 
1 .9938. The resulting quantized, scaled symbols have 5 bits for each component. 

It should be appreciated that any of the foregoing methods may be tangibly 
embodied as a series of instructions stored on a processor readable medium or 
memory including, without limitation, RAM, ROM, PROM, EPROM, EEPROM, 
hard disk, floppy disk, CD-ROM, DVD, tape, flash memory, etc. A system may also 
be formed comprising the processor readable medium coupled to a processor 
configured to access and execute the instructions stored on the processor readable 
medium. 

A system for performing the method may also be designed and executed using 
synthesized (digital) combinatorial logic elements, arithmetic elements, and memory. 
System control may be implemented by a finite state machine embodied in the form 
of ASIC hardware. The ASIC hardware may be synthesized using combinatorial 
logic, e.g. AND/OR/NAND/NOR/XOR/NOT gates and flip-flops. In this 
implementation, the finite state machine may be used for control, RAM may be used 



SCALE = 10^'''''''-'^' 



(15) 
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for table storage, and arithmetic xmits (adders, multipliers) may be used to perform the 
necessary arithmetic for interpolation, scaling, etc. In one implementation example, 
software chip design tools may be used to synthesize the algorithm as digital logic, 
with state machine control. In this implementation example, there may be no distinct 
5 CPU on the actual chip that formally processes 'instructions'. However, it should be 
appreciated that a system in which a CPU formally processes instructions is still 
within the scope of the invention. 



Any of the methods of the invention may be embodied as hardware, software, 
10 or a combination of hardware and software. The hardware may comprise dedicated 
analog or digital circuitry, integrated circuits, ASICs, PLAs, or the like. The software 
may comprise instructions executable by a DSP. For purposes of this disclosure, the 
term "logic" will encompass hardware, software, or a combination or hardware and 
software for performing a function. 

15 While various embodiments of the invention have been described, it will be 

apparent to those of ordinary skill in the art that many more embodiments and 
implementations are possible that are within the scope of this invention. Accordingly, 
the invention is not to be restricted except in light of the attached claims and their 
equivalents. 
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